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O . Abstract 

o ■ 

""^ ■ We classify all spherically symmetric dust solutions of Einstein's equa- 

Q H . tions which are self-similar in the sense that all dimensionless variables 

$_i \ depend only upon z = r/t. We show that the equations can be reduced 

jpf)" to a special case of the general perfect fluid models with equation of state 

J> ■ p = aa. The most general dust solution can be written down explicitly 

^ , and is described by two parameters. The first one (E) corresponds to the 

^ \ asymptotic energy at large \z\, while the second one (D) specifies the value 

of z at the singularity which characterizes such models. The E = D = 
solution is just the flat Friedmann model. The 1-parameter family of so- 
lutions with z > and D = are inhomogeneous cosmological models 
which expand from a Big Bang singularity at t = and are asymptoti- 
cally Friedmann at large z; models with E > are everywhere underdense 
relative to Friedmann and expand forever, while those with E < are 
everywhere overdense and recollapse to a black hole containing another 
singularity. The black hole always has an apparent horizon but need not 
have an event horizon. The D = solutions with z < are just the time 
reverse of the z > ones. The 2-parameter solutions with D > again 
represent inhomogeneous cosmological models but the Big Bang singular- 
ity is at z = —1/D, the Big Crunch singularity is at z = +1/D, and any 
particular solution necessarily spans both z < and z > 0. While there 
is no static model in the dust case, all these solutions are asymptotically 
"quasi-static" at large \z\. As in the D = case, the ones with E > 
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expand or contract monotonically but the latter may now contain a naked 
singularity. The ones with E < expand from or recollapse to a second 
singularity, the latter containing a black hole. The D < models either 
collapse to a shell-crossing singularity and become unphysical or expand 
from such a state. 

PACS numbers: 0420, 0420J, 9530S, 9880H 

1 Introduction 

Spherically symmetric self-similar solutions to Einstein's equations have the fea- 
ture that every dimensionless variable is a function of some dimensionless combi- 
nation of the cosmic time coordinate t and the comoving radial coordinate r. In 
the simplest situation, a similarity solution is invariant under the transformation 
r — > at,t — > at for any constant a, so the similarity variable is z = r/t. Geo- 
metrically this corresponds to the existence of a homothetic Killing vector. Such 
solutions have been the focus of much attention in General Relativity because the 
field equations simplify to ordinary differential equations (Cahill & Taub 1971). 

Even greater simplification is afforded if one focusses on the situation in which 
the source of the gravitational field is pressureless dust since, in this case, the 
solutions can often be expressed analytically and are just a special subclass of 
the more general spherically symmetric Tolman-Bondi solutions (Tolman 1934, 
Bondi 1947, Bonnor 1956). A number of people have studied such solutions: 
for example, Gurovich (1967), Carr & Hawking (1974), Dyer (1979), Wu (1981), 
Maharaj (1988), Carr & Yahil (1990), Ori & Piran (1990), Henriksen & Patel 
(1991) and Sintes (1996). These papers are described in some detail by Krasinski 
(1997). As examples of applications of particular physical intererst, self-similar 
dust solutions have also been used in studying naked singularities by Eardley & 
Smarr (1984), Christodoulou (1984), Lake (1992), Joshi & Dwivedi (1993), Joshi 
& Singh (1995), Dwivedi & Joshi (1997) and in modelling cosmic voids by Tomita 
(1995, 1997a, 1997b) and Carr & Whinnett (1999). 
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In this paper we will present a classification of spherically symmetric homo- 
thetic dust models which is "complete" subject to certain specified restrictions. 
This serves as a first step in the more general analysis, presented in an accompa- 
nying paper (Carr & Coley 1999a), of spherically symmetric self-similar perfect 
fluid models with equation of state p = a/i. Despite the simplifications entailed 
in dropping pressure, we will find that many of the features of the dust (a = 0) 
solutions carry over to the q^O case, at least in the supersonic regime. In par- 
ticular, many (though not all) of the types of solutions with pressure have direct 
analogues in the dust case. Our aim is therefore to use the exact analytic dust so- 
lutions to derive qualitative features of self-similar solutions which will also turn 
out to pertain when there is pressure but which can then only be demonstrated 
numerically. It should, of course, be stressed that the introduction of pressure 
is also associated with many new features - especially in the subsonic regime - 
which cannot be understood in this way. 

The plan of the paper is as follows: In Section 2 we will show how the dust 
equations can be regarded as a special case of the equations with pressure provid- 
ing one adopts a specific prescription in taking the limit a — > 0. In Section 3 we 
discuss the general family of self-similar solutions, showing that they can be con- 
veniently categorized as asymptotically Friedmann and what we term asymptot- 
ically quasi-static. In each case, we will focus on the important physical features 
of the solutions, such as the presence of apparent horizons, event horizons and 
singularities. We will also clarify the connection between models with positive 
and negative z. We will specify the sense in which our classification is "complete" 
in Section 4. 

2 Spherically Symmetric Similarity Solutions 

In the spherically symmetric situation one can introduce a time coordinate t 
such that surfaces of constant t are orthogonal to fluid flow lines and comoving 
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coordinates (r, 0, <fi) which are constant along each flow line. The metric can then 
be written in the form 



ds z = e 2u dt 2 - e 2X dr 2 - R z dtf, dQ' z = d0 l + sin^ , 



(2.1) 



where u, A and R are functions of r and t. The Einstein equations have a first 
integral 



m r 



\R 



-2A 



<9t / y <9r 

This can be interpreted as the mass within comoving radius r at time t: 

,dR 



m ( r ) = 4n / uR — — dr . 
Jo or 



(2.2) 



(2.3) 



where /x(r, t) is the energy density and we choose units in which c = G = 1. This 
is constant if p — 0. Eqn (2.2) can be written as an equation for the energy per 
unit mass of the shell with comoving coordinate r: 

*.^-D-^-=. —(f), r^(f). (,4) 

This can be interpreted as the sum of the kinetic and potential energies per unit 
mass. E and V are conserved along fluid flow lines in the p = case. 

By a spherically symmetric similarity solution we shall mean one in which the 
spacetime admits a homothetic Killing vector £ that satisfies 



(2.5) 



This means that the solution is unchanged by a transformation of the form t — > at, 
r — > ar for any constant a. Solutions of this sort were first investigated by Cahill 
& Taub (1971), who showed that by a suitable coordinate transformation they 
can be put into a form in which all dimensionless quantities such as u, A, E and 



R 



M 



m 
R/ 



P = pR 2 , W = \iR 2 



(2.6) 
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are functions only of the dimensionless variable z = r/t. This means that the 
field equations reduce to a set of ordinary differential equations in z. Another 
important quantity is the function 



which represents the velocity of the surfaces of constant z relative to the fluid. 
These surfaces have the equation r = zt and therefore represent a family of 
spheres moving through the fluid. The spheres contract relative to the fluid for 
z < and expand for z > 0. This is to be distinguished from the velocity of the 
spheres of constant R relative to the fluid: 



This is positive if the fluid is collapsing and negative if it is expanding. Special 
significance is attached to values of z for which |V| = 1 and \Vr\ = 1. The 
first corresponds to a Cauchy horizon (either a black hole event horizon or a 
cosmological particle horizon) and the second to a black hole or cosmological 
apparent horizon. One can show that the existence of an apparent horizon is also 
equivalent to the condition M = 1/2. 

Although our main focus in this paper is (pressureless) dust solutions, it is 
elucidating to start by considering the equations for a fluid with pressure. The 
only barotropic equation of state compatible with the similarity ansatz is one of 
the form p = afi (—1 < a < 1). If one introduces a dimensionless function x(z) 
defined by 



V(z) 




(2.7) 




(2.8) 



x 



(z) = {Anfir 2 )- a/il+a \ 



(2.9) 



then the conservation equations T^ u . u = can be integrated to give 



e u = (3xz 



2a/(l+a) 



(2.10) 




1/a c2 



(2.11) 



5 



where (3 and 7 are integration constants. The remaining field equations reduce 
to a set of ordinary differential equations in x and S: 

S + S+ f-|-| - lf\ [S + (1 + a)S] = 0, (2.12) 
yl + ao ax J 

( 2a ^ 2 ^ <?4 + 2 5 (2 _ 2Q)/az(2 _ 2Q)/(1+Q) _ 7 ^ 4 i _ x \ = (1 + a)a .(i-« ) /« j 



1 + a J (3 2 S x \ a 

M = S 2 x- { - 1+a)/a 



1 + + 



(2.13) 
(2.14) 



M 2 



1 + ^Lx-V^)/^)^ 2 - ^yV^S 8 (l + I) , (2.15) 
where the velocity function is given by 

V = ^-i x (i-^)/a S -2 z (i-a)/(i + a) (2.16) 

and an overdot denotes zd/dz. The other velocity function is 

VS 

V R =— ^ (2.17) 

s + s v ; 

and, from eqns (2.4) and (2.15), the energy function is 

E = \l 2 *- {2,a) S & (l + S 1 ^ -\ , (2.18) 

which necessarily exceeds —1/2. 

As discussed by Carr & Yahil (1990), we can best envisage how these equations 
generate solutions by working in the 3-dimensional (x, S, S) space. At any point 
in this space, for a fixed value of a, eqns (2.14) and (2.15) give the value of 
z; eqn (2.13) then gives the value of x unless \V\ = ^fa and eqn (2.12) gives 
the value of S. Thus the equations generate a vector field (x, S, S) and this 
specifies an integral curve at each point of the 3-dimensional space. Each curve is 
parametrized by z and represents one particular similarity solution. This shows 
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that, for a given equation of state parameter a, there is a 2-parameter family 
of spherically symmetric similarity solutions. In general there would be a sonic 
point, with possible associated discontinuities, at \V\ = fa. This corresponds 
to crossing a 2-dimensional surface in the solution space but we need not discuss 
this complication here. 

We now focus on the dust solutions. Although eqns (2.10) to (2.16) break 
down when a = 0, in that some of the terms appearing there disappear or diverge, 
we will show that they are still formally applicable providing the function x 
defined by eqn (2.9) is set to 1 whenever it does not appear with the exponent 
1/a. Otherwise one must make the substitution 

xV- _ (4^r 2 )- 1 , (2 .i9) 
a x dlnz 

as suggested by eqn (2.9). In fact, the relevant equations are most simply obtained 

by noting that both the energy and mass within comoving radius r are conserved, 

so that E and m/r = MS are constant. If we put m/r = k, then eqn (2.3) implies 

^ndR dm _ ^ . 

AttuR 2 — = — = k 2.20 
or dr 

and this can be combined with eqns (2.9) and (2.11) to give 

e A = 7" 1 (47r/ir 2 5 2 )- 1 = k' 1 ^ 1 ^. (2.21) 

or 

On the other hand, eqn (2.4) implies e A = Y~ l (dR/dr) where V = ±^1 + 2E 
and so the constant k is just r/7. The mass function is therefore 

where we have taken the positive square root for T to ensure that the mass if 
positive. (We discuss the negative mass case later.) Putting x — 1 in eqn (2.10) 
also gives e v = (3, so eqn (2.4) can be rewritten as 



dS 

dz 



2 



VI + 2E 

is ■ 



(2.23) 
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Once this equation has been integrated to give S(z), all the other functions can 
be obtained, so one has solved the problem completely. 

We now show that eqns (2.10) to (2.16) are all formally satisfied if one uses 
the prescription given by eqn (2.19). Eqn (2.14) can be written as 

,dR 



MS = m/R = S 2 (47ri2r 2 )(S + S) = Ait/jR 2 



dr 



and this is just equivalent to eqn (2.20). Since eqn (2.21) implies 

(47r/ir 2 )- 1 = 1 S 2 (S + S)/T, 



we also have 



dln{jir 2 ) 



2S_ _ S + S 



(2.24) 



(2.25) 



(2.26) 



dlnz S S + S 

Eqn (2.12) is then automatically satisfied and from eqn (2.23) formally corre- 
sponds to E — 0. Finally we can substitute for the x terms in eqn (2.13), using 
eqns (2.19) and (2.25), to obtain 

r 



s + s 



S 2 z 2 



(2.27) 



and this can be integrated to give eqn (2.23). 

It is now convenient to scale the r and t coordinates so that (5 — 7 = 1. Eqn 
(2.23) then implies 

dS_ ^2E + 2T/S 



dz 

and this can be integrated to give 



(2.28) 



z 



VES 2 +rs 

V2E 



2F 



^sinh-Vf (E>0) 



-1 l -ES 4. \ZES' 2 +TS 

V T y/2E 



(E = 0) 

(-1/2 < E < 0) 



(2.29) 



L (-2£)3/ 2 S111 

where D is an integration constant and sin -1 is taken to lie between and tt. In 
the first two cases, the upper and lower signs on the left apply for dS/ dz positive 
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and negative, respectively, this sign being constant for any particular solution. In 
the third case, the sign on the left is fixed for a particular solution, even though 
dS/dz may switch sign, but the sign of the last term is plus if sin -1 lies between 
and 7r/2 and minus if it lies between 7r/2 and ir. If we took the negative square 
root in eqn (2.22), corresponding to M and T being negative, there would be 
another solution for E > given by 

1 Jes 2 - \T\S 2|r| , [ES 

D ^- z - V2E + (d^ CQSh Vw {E>0) - (2 - 30) 

This solution is unphysical, since the mass is negative, but it is of interest for 
comparison with the solutions with pressure. 

Eqns (2.16), (2.17), (2.19), (2.25) and (2.28) give the velocity functions as 



Sz ± J2E + 2T/S 
V= V— (2.31) 



and 



'2E + 2T/S 

V R = ±* f , (2.32) 



while eqns (2.16) and (2.31) imply that the density is given by 

4^t 2 = -4^7 = F (2.33) 

Z S 2 V zS i(Sz ± ^2E + 2T/S) 

where the upper and lower signs again correspond to dS/dz being positive and 
negative, respectively. Note that V and ji is negative (corresponding to tachyonic 
models) for the solution given by eqn (2.30) and \x is also negative. In all cases 
the metric can be written as 

1 + 2E ' v ; 

which is the standard Tolman-Bondi form with constant energy function E(r). 
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3 Classification of Solutions 



Eqn (2.29) implies that there is a 2-parameter family of similarity solutions (as 
in the general a case). In this section we will provide a complete description of 
these solutions. We will start by considering the simplest one: the flat Friedmann 
solution with D — E — 0. We will then consider the one-parameter family of 
solutions with E ^ 0, D = 0. Finally we will consider the full two-parameter 
family of solutions with £ ^ 0,D / 0. In each case, we will show the form of 
the physically interesting quantities S, V and fit 2 as functions of z. In obtaining 
the full family of solutions, it is crucial that we allow z to be either positive or 
negative. Our analysis will also cover the (presumably unphysical) solutions with 
negative mass because they relate to some of the solutions with pressure. 

3.1 E = D = solution 

In this case eqns (2.21), (2.29), (2.31) and (2.33) give 

S = (V2z/3)' 2/3 = M~\ V = (z/Q) 1/3 , fi = (Qui 2 )- 1 . (3.1) 

This corresponds to the standard dust Friedmann model with zero curvature 
constant. The metric can be put in the usual form by making the substitution 
f = (9r/2) 1 / 3 , which gives 

ds 2 = dt 2 -t 4 ^[dr 2 + f 2 dn 2 ]. (3.2) 

Note that the curvature constant must be zero because otherwise there would be 
an intrinsic scale, which would contradict the similarity assumption. 

3.2 D = solutions 

Solutions with D = are asymptotically Friedmann as \z\ — > oo and are specified 
entirely by the energy parameter E and were studied by Carr & Hawking (1974) 
and Carr & Yahil (1990). The form of S(z) in these solutions is shown in Figure 
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(la), the arrows always corresponding to the direction of increasing time. The 
solutions with z > correspond to initially expanding Big Bang models: they 
start from an initial Big Bang singularity (S — 0) at t — [z — oo) and then 
either expand indefinitely (S — > oo) as t — > oo (z — > 0) for E > or recollapse to 
a black hole singularity (S — 0) at 

(-2gff 

Z5 = WTTH (3 - 3) 

for < 0. Note that z$ corresponds to the physical origin since R = rS = 
there. The form of V(z) in the z > solutions is shown in Figure (lb). In the 
first case, V decreases monotonically from oo to 0. In the second case, it reaches 
a minimum before rising to oo at zs- One can show that the values of z and V 
at the minimum both decrease as E increases; the minimum will exceed 1 (in 
which case the whole Universe is inside the black hole) if E is less than some 
critical negative value E* and it will be less than 1 (in which case there is a black 
hole event horizon and a cosmological particle horizon) if E exceeds E*. The 
solutions with z < are the time-reverse of the z > ones and the sign of V is 
also reversed: as t increases from — oo to (i.e. as z decreases from to — oo), 
the E > models collapse from an infinitely dispersed state (S = oo) to a Big 
Crunch singularity (S = 0); the E < models also collapse to a Big Crunch 
singularity but they emerge from a white hole and are never infinitely dispersed. 

Both S and V have the same z-dependence as in the E = Friedmann solution 
as \z\ — > oo: 

S w [9VTT2E/2} 1 ^\z\~ 2 / 3 , V w [6(1 + 2E)]- 1 / 3 z 1 ^. (3.4) 

However, the E ^ solutions deviate from the E — solution at small values of 
\z\. The E < solutions never reach z = at all, while the E > ones have 

S ^ {2Ef/ 2 \z\-\ V ^-{l + 2Ef/ 2 E- 1 z\n[{2Ef' 2 (l + 2E) 1 / 2 \z\} (3.5) 
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as \z\ — > 0. The first relation implies that the circumference function R(r, t) = Sr 
is non-zero in limit r — > unless = since 

R(0, t) = V2E t. (3.6) 

This means that the "coordinate" origin (r = 0) is an expanding 2-sphere. [This 
feature is specific to the dust case and does not arise if there is pressure.] This 
has a natural physical interpretation since the forms of S and V are similar to 
those in the Kantowski-Sachs (1966) solution, in which all the matter is localized 
on a shell [cf. eqns (3.19) and (3.23) in Carr & Coley (1999a)]. However, we 
note that there is no exact self-similar Kantowski-Sachs solution in the dust case. 
To obtain a complete solution, one must therefore match the self-similar solution 
onto a (non-self-similar) part inside R(t, 0). In the E < case, we have seen that 
the physical origin is the black hole singularity zs, so only for E = can one 
identify z = with the physical origin. 

The form of the density function fit 2 can be derived from eqn (2.33) and is 
shown in Figure (lc). For a given fluid element, this specifies the density as 
a function of time fi{t). For a given time, it also specifies the density profile 
/i(r) and this illustrates that a non-zero value of E necessarily introduces an 
inhomogeneity into the model. Solutions with E > are everywhere underdense 
relative to the Friedmann model, with eqns (2.33) and (3.5) implying that fit 2 
goes to as (In l^l) -1 as z — > 0. (This suggests that the interior non-self-similar 
region should be a vacuum.) Solutions with E < are everywhere overdense 
relative to Friedmann, with fi diverging at the singularity. Note that eqns (2.33) 
and (3.4) imply that fit 2 is independent of E to 1st order as \z\ — > oo. 

The form of the mass function M(z) in the D = solutions is not shown 
explicitly but can be immediately deduced from the expression for S since eqn 
(2.21) gives M = T/S. In the E > case, there is always a single point where 
M — 1/2 and this corresponds to the cosmological apparent horizon. In the 
E < case, eqn (2.29) implies that S has a maximum of r/]i?| and so eqn (2.21) 
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shows that M has a minimum of \E\. Since this is less than 1/2, there are always 
two points where M — 1/2, one corresponding to the black hole apparent horizon 
and the other to the cosmological apparent horizon. Note that a black hole's 
apparent horizon always lies within or coincides with its event horizon (Hawking 
& Ellis 1973), which is why the first can exist without the second. Eqns (2.21) 
and (3.3) imply that the mass associated with this singularity is 

m s = (MSz)s t = (-2£) 3/2 t/(2?r). (3.7) 

It therefore starts off zero when the singularity first forms at t — but then grows 
as t. The mass of the black hole is given by a similar formula but with z having 
the value appropriate for the event horizon or apparent horizon. Since the former 
may not exist, it is more appropriate to use the latter. 

Finally, we consider the E > negative-mass solutions given by eqn (2.30). 
Their form is indicated by the dotted curves in Figures (1). S, V and // have the 
same form as in the positive mass solutions for small values of \z\ except that V 
and \x reverse their signs. However, the solutions are very different at large values 
of \z\ since eqn (2.30) shows that S must always exceed \F\/E. Indeed it tends 
to this value asymptotically, so we have 

S « VI + 2E/E, V « —z/E, fir 2 « -E 3 /y/l + 2E (3.8) 

as \z\ — > oo. The form of this solution is closely related to that of the a « 1 
static solution [cf. eqn (3.29) of Carr & Coley (1999a)], although there is no static 
solution in the a = case itself. 

3.3 E = solutions 

We now put E — and consider the effect of introducing a non-zero value for 
the constant D. [In this case, eqn (2.28) does not permit T < 0, so there are no 
negative-mass solutions.] The form of S(z) for the D > solutions is shown in 
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Figure (2a). There are two types of solutions in this case, one expanding and the 
other collapsing. For the expanding solutions (solid lines), S = aX z = —1/D 
and so the Big Bang occurs before t = (i.e. it is "advanced"). As t increases to 
(i.e. as z decreases to — oo), S tends to the finite value 

Soo(D) = (3D/V2) 2 / 3 . (3.9) 

As t further increases from to +oo (i.e. as z jumps to +oo and then decreases 
to 0), S increases monotonically to oo. For the contracting solutions (broken 
lines), S starts infinite at t — — oo [z = 0) and then decreases to S^D) as t 
increases to [z — > — oo). As t further increases (i.e. as z jumps to +oo and 
then decreases), S continues to decrease until it reaches at the Big Crunch 
singularity at z — 1/D. Both types of solutions are characterized by the fact that 
they have just one singularity and span both positive and negative values of z. 
Note that for each value of the two asymptotic solutions just correspond to 
the plus and minus signs in eqn (2.28). 

The form of V(z) in the D > solutions is shown in Figure (2b). For the 
expanding solutions (solid lines), it starts off at — oo at the Big Bang (z = —1/D), 
reaches a negative maximum and then, from eqn (2.31), tends to 

V = S OQ (D)z = (3D/V2) 2/3 z (3.10) 

as z — > — oo. When z jumps +oo, V becomes positive but eqn (3.10) still applies. 
As z decreases from +oo to 0, V decreases monotonically to 0. For the contracting 
models (broken lines), V starts from zero at z = and monotonically decreases as 
z goes to — oo, being again given by eqn (3.10) asymptotically. When z jumps to 
+oo, V jumps to +oo and then decreases to a minimum before rising to infinity at 
the Big Crunch singularity. Note that the maximum value of V for the expanding 
solutions will exceed —1 and the minimum value for the contracting ones will be 
less than +1 (i.e. the minimum value of |V| is less than 1) if D exceeds some 
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critical value D + . A simple calculation (see later) shows that the values of \z\ 
and |V] at the stationary point are given by 

- (2+ ^ iH„,„4^±^y /3 , (3,D 



\^\min I r^j-y J i | r \min I I 

so the stationary point moves towards the origin as D increases and D + = 26/3 + 
5\^3 ~ 17. For D > D + , the condition \V\ = 1 will be satisfied at three values 
of z. As illustrated by Figure 14 of Ori & Piran (1990), this means that the 
contracting solutions will form a black hole in which the central singularity is 
naked for a while. 

The crucial feature of these solutions is that, while the form of V(z) is like 
that in the D = case for small \z\, V scales as z rather than z 1 ^ [cf. eqn 
(3.4)] for large \z\. This is because any solution with finite S at infinity must be 
"nearly" static in the sense that dS/dz tends to zero. However, the solutions are 
not asymptotic to an exact static solution (indeed this does not exist in the a = 
case) because eqn (2.32) implies that Vr tends to a non-zero value: 

^ = ±(^) • (3-12) 

We therefore term these solutions asymptotically "quasi-static" . If V£° is positive, 
the fluid is collapsing at infinity; if VJf is negative, it is expanding. Note that eqn 
(2.28) implies that both dS/dz and zdS/dz tend to zero at large \z\ but z 2 dS/dz 
[which directly relates to from eqns (2.28) and (2.32)] tends to a non-zero 
value except in the limit D — > oo. 

The form of the density function [it 1 in the D > solutions is also interesting 
and is illustrated in Figure (3c). From eqn (2.33) the density parameter is given 
by 

" B ^ = (1 T 3^1(1 (3 - 13) 
where the upper and lower signs apply for positive and negative values of dS/ dz, 
respectively. [The inverse of the factor 6Trt 2 corresponds to the density in a flat 
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Friedmann dust universe, as indicated by eqn (3.1).] For a given fluid element, 
this describes how the density evolves as a function of time and it has the expected 
form. However, at a given time it also prescribes the density profile and one sees 
immediately that a non-zero value of D (like a non-zero value of E) introduces 
an inhomogeneity. This inhomogeneity has a particularly interesting form. In 
the z < regime, the profile for the collapsing solutions is homogeneous for 
\z\ « 1/D but has /i ~ r~ 2 for \z\ >> 1/D (i.e. it resembles an isothermal 
sphere with a uniform core). In the z > regime, the collapsing solutions again 
have /i ~ r -2 for \z\ » 1/D but the density diverges at z — 1/D (i.e. one has 
a density singularity at the centre of an isothermal sphere). For the expanding 
solutions, the signs of z are reversed. These features are illustrated in Figure (2c) 
and have an obvious physical interpretation. 

Although the asymptotically quasi-static solutions have an natural cosmolog- 
ical interpretation when z is allowed to span both positive and negative values, 
we see that the z > and z < solutions also have a non-cosmological interpre- 
tation when considered separately: they just represent collapsing and expanding 
self-similar models which evolve from an initially isothermal distribution. It is 
interesting that the isothermal model (which is usually associated with a static 
solution) features prominently in both regimes, despite the fact that there is no 
exact static solution in the dust case. Note that the behaviour at z = is different 
from the D = case, in that Q tends to 1 rather than for E > 0. 

The mass function is M = S^ 1 in this case and therefore decreases or increases 
monotonically. There is just one value of z at which M — 1/2 (corresponding 
to a black hole or cosmological apparent horizon) but this may be in either the 
positive or negative z region. Eqn (3.9) implies that the asymptotic value of M as 
\z\ — > oo is less than 1/2 for D > 4/3. In this case, the collapsing solutions have 
their apparent horizon in z > 0, whereas the expanding ones have it in z < 0. 
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The mass of the (possibly naked) singularity in these solutions is 

m s = (MSz)s t = t/D (3.14) 

from eqn (2.21). As in the asymptotically Friedmann case, it starts off zero at 
t — but then grows as t. In the limit D — > oo, one gets a naked singularity of 
zero mass at the origin (cf. the static solution). 

Finally, we consider the D < solutions. The form of S(z) in this case 
is shown by the dotted curves in Figure (2a). Such solutions are confined to 
\z\ < —l/D, with S either decreasing monotonically for z < (i.e. as t increases 
from — oo) or increasing monotonically for z > (i.e. as t increases to +oo). 
However, these solutions break down when S is too small. This is because eqn 
(2.31) implies that |V| increases to some maximum value and then falls to zero 
at \z\ = —1/(3D); this is indicated by the dotted curve in Figure (2b). From eqn 
(2.33) this means that the density diverges there, as shown by the dotted curve in 
Figure (2c). This divergence is associated with the formation of a shell-crossing 
singularity since the model resembles the Kantowski-Sachs solution at this point. 
For —l/D > \z\ > —1/(3D), the density and velocity functions become negative 
but this is presumably unphysical. 

3.4 D^0,£/0 solutions 

The forms of S(z) for the (D > 0, E ^ 0) solutions are indicated in Figure (3a). 
The figure assumes that D is fixed but allows E to vary. The (D > 0, E > 0) 
solutions are qualitatively similar to the (D > 0, E — 0) ones in that they are 
monotonically expanding or collapsing and span both positive and negative z, 
as illustrated by the upper solid and broken curves, respectively. They are also 
asymptotically quasi-static, in the sense that S tends to a finite value as \z\ — > oo, 
even though Vr is non-zero there from eqn (2.32). The form of the solutions near 
z = is still given by eqns (3.5) for E > 0, so the behaviour is like that in 
the (E > 0,D = 0) case here. In particular, z = no longer corresponds to 
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the physical origin, so one again has to attach the solution to a non-self-similar 
central region. 

The (D > 0, E < 0) solutions are qualitatively different from the (D > 0, E = 
0) ones in that the models no longer collapse from or expand to infinity. This is 
clear from eqn (2.28), which implies that S has a maximum value of T/\E\, so 
all the solutions start off expanding and then recollapse. Note that there is no 
exact static solution since that would be incompatible with eqn (2.13), the term 
on the right-hand-side being non-zero for a = 0. Figure (3a) shows that there are 
two types of (D > 0, E < 0) solutions. One type (illustrated by the lower solid 
curves) expands from the Big Bang singularity at z — — 1/D and then recollapses 
to a black hole singularity at 

-i -i 



zs 



2nVT+2E 



(3.15) 



(-2E) 3 / 2 

This reduces to the value given by eqn (3.3) if D = 0. The other type (illustrated 
by the lower broken curves) expands from a white hole singularity at — zs and 
then recollapses to a Big Crunch singularity at z — 1/D. 

Eqn (3.15) implies that zs = 1/D, so that the solution is symmetric in z, if 
D has the value 

^sym- (_2£)3/2 " [6.10) 

One can invert this condition to obtain the associated value of E in terms of D: 

(3.17) 



47T 

E svm = = — sinh 

y V3D 



1 , _i {3V3D\ 
- sinh ; 



V 8n J. 

and this specifies a 1-parameter family of solutions with = 0. If one considers 
the limit of the symmetric solution as D — > 0, one finds E sym — > —1/2 and 
zs — > oo. On the other hand, if one considers the limit as D — > oo, one finds 
E S ym and zs — > 0, so that both singularities go the origin. This is the closest 
one can get to a static solution in the dust case. Note also that zs — > in the 
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limit E — > 0~ whatever the value of D; the sudden transition as one goes from 
E = 0~ to E = + is illustrated in Figure (3a). 

This value of E given by eqn (3.17) has a special physical significance in that 
it prescribes the minimum value of E allowed for given D, as indicated by the 
lower boundary in Figure (4). The proof of this is as follows. If one takes the 
limit of eqns (2.14) and (2.15) with a = as z — > oo, using eqns (2.19) and (2.32) 
and the fact that S — > from eqn (2.28), one obtains 

z 2 S* = (1 + 2E){V™? = (4^r 2 )^ + (Stt^ 2 )^ - 1 (3.18) 

where one can regard (fir 2 )^ and as independent asymptotic parameters. If 
one now fixes (/xr 2 )^ and assumes that it is positive, then the right-hand-side of 
eqn (3.18) decreases monotonically with decreasing S^- One therefore gets a real 
solution for Vr only if Soo exceeds a certain value and - by monotonicity - this 
must be the value associated with the symmetric solution. Thus a real solution 
(with positive density) exists only for E > E sym or, equivalently, D < D sym . This 
means that z$ is always positive and less than 1/D and that the maximum of S 
will always occur at the opposite sign of z as the \z\ = 1/D singularity. These 
features are illustrated by the curves in Figure (3a). 

The form of V(z) in the these solutions is shown in Figure (3b). For (D > 
0, E > 0) it is similar to that in the (D > 0, E — 0) case. As z decreases from 
to — oo, V decreases monotonically from to — oo; it then jumps to z — +oo 
and, as z continues to decrease, it falls to a minimum and rises to infinity at 
the Big Crunch singularity at z — 1/D. Also as in the (D > 0, E = 0) case, 
this minimum will fall below 1, corresponding to a naked singularity, providing 
D exceeds some value D + (E). We derive an implicit expression for D + (E) later 
but, for the present, note that it increases with increasing E and reduces to the 
value D + which arose in Section 3.3 when E — 0. The condition for a naked 
singularity can also be expressed as the requirement that E exceed some critical 
value E+(D). 
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For (D > 0, E < 0) the form of V is similar to that in the (D = 0, E < 0) 
case. As z decreases from —1/D, V rises from — oo until some maximum and then 
falls to — oo quasi-statically as z — > — oo. It then jumps to z — +oo and falls to a 
minimum before rising to oo at z s . As in the (D — 0, E > 0) case, there will be 
a black hole event horizon if the minimum value of V is less than 1; this requires 
that E exceed some value E*(D), which must reduce to the value E* given in 
Section 3.2 when D — 0. We derive an implicit expression for E*(D) below. It 
should be stressed that the value E*(D) is associated with the minimum of |V| 
near the singularity at \z\ = zs and is different from the value E + (D) associated 
with the minimum near \z\ = 1/D. The relationship between the values E + (D) 
and E*(D) is discussed below. 

In order to understand the form of the curves in Figures (3a) and (3b) more 
precisely, it is useful to specify their asymptotic behaviour. As \z\ — > oo, eqns 
(2.29) and (2.31) with E > imply that S(z) and V(z) have the following asymp- 
totic forms: 

/ 2E \ 

» DV2E, \4o w {y+2EJ DZ (3,19) 
for D » [(1 + 2E)/E 3 } 1 / 2 and 

for D « \{1 + 2E)/ E z ] 1 / 2 . The transition value for D between these two regimes 
is just an extrapolation of the expression for D sym given by eqn (3.16) into the 
E > regime; it scales as E' 1 for E » 1 and E~ 3/2 for E « 1. Note that 
eqn (3.20) agrees with eqns (3.9) and (3.10) in the limit E = 0. For £ < 0, eqn 
(3.20) still applies if D « D sym but one has 

S OQ ^Vl + 2E/\E\, V^^z/lEl (3.21) 

for D ?a -D sym (i.e. tends to the value associated with the symmetric solution). 

These equations prescribe the asymptotic forms for an d in the different 
(E, D) regimes of Figure 4. For fixed D, eqns (3.19) to (3.21) show that always 
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increases with E but is roughly constant for \E\ « 1. The behaviour of Voo is 
more complicated: for D « 1, it first decreases with increasing E and then 
flattens off; for D » 1, it first increases with increasing E before flattening off. 
These features are indicated in the Figures (3a) and (3b). Note that all these 
solutions are quasi-static and not exactly static asymptotically since eqn (2.32) 
gives 

4 1 / 3 / 2\E\ \ 1/2 

lT ^ ± (3D)V3(i + 225)V3' V k =± [TT2e) ^ 

for D » [(1 + 2E)/E 3 Y/ 2 and D « [(1 + 2E)/\E\ 3 } 1 / 2 . The sign is positive 
for collapsing solutions and negative for expanding ones. Note that the first 
expression agrees with eqn (3.12) in the limit E — 0. 

We now derive implicit expressions for the functions E*(D) and E + (D), which 
are related to the existence of event horizons or naked singularities. Differentiat- 
ing eqn (2.31) shows that when dV/dz = one always has 

V = l/{S 2 z) (3.23) 

and eqn (2.31) then gives 

VT - 1/(VS) = ±^2E + 2Y/S , (3.24) 

where the positive and minus signs corresponds to the sign of dS/dz. If one also 
requires |V| = 1 at the stationary point, eqn (3.24) gives two roots 

s = 2T ± V4r 2 - 1, |*| = sr 2 - 1 ^4rV4r 2 - 1, (3.25) 

where the plus and minus signs are distinct from the ones appearing in eqn (3.24). 
These roots can be real providing T > 1/2, corresponding to E > —3/8. However, 
one needs to check whether both of these solutions satisfy condition (3.24). 

Inserting the solutions (3.25) into eqn (2.29) gives an expression for D in 
terms of E. This expression is complicated in general but it simplifies in certain 
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regimes. For E » 1, which implies V 



\plE >> 1, one obtains two possible 



solutions: 



S^AT, \z\t* 



i/(ier 2 ), 



D as 32E 



(3.26) 



and 



5wl/(4r), 




lnv^l « 0.06E 



-i 



(3.27) 



The first has T > 1/S and therefore requires dS/dz > from eqn (3.24), which 
leads to a consistent solution. However, the second has T < 1/S and requires 
dS/dz < 0, which does not. In the limit E — > 0, which implies fwl, one obtains 



The upper sign gives T > 1/S and therefore requires dS/dz > 0, which leads to a 
consistent solution as E tends to from either above or below. Eqn (2.29) then 
gives 



Note that the constant part of this expression gives the same limiting value of 
D as implied by eqn (3.11). The lower sign in eqn (3.28) gives T < 1/S and 
requires dS/dz < 0, which does not lead to a consistent solution. Finally we 
note that eqns (3.25) and (2.29) lead to a unique value of E and D for which the 
symmetric solution has V m in = 1; this necessarily corresponds to a point on the 
lower boundary in Figure 4. 

These limiting behaviours allow one to infer the rough form of the functions 
E*(D) and E+(D), as indicated in Figure (4). Here we have used the fact that 
dE + /dD is positive as E — > 0, as follows from eqn (3.29). The form of the 
functions in the E < regime can be inferred from the fact that E*(D) must 
reach the value E* mentioned in Section 3.2 when D = (although this value 
has not been calculated explicitly). Also E + (D) and E + (D) must reach the line 
E = E sym (D) at the same value of E and this must clearly exceed —3/8 from 



S w (2±V3)(1±2E/V3), 



z ~ 



(7^4^(1^4^/73). 



(3.28) 




(3.29) 
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eqn (3.25). We note that, for sufficiently large values of D, there may be both an 
event horizon and a naked singularity. 

The form of /it 2 in these solutions is shown in Figure (3c), although this gives 
only some of the solutions shown in Figures (3a) and (3b). It can be understood 
as a composite of the curves shown in Figure (lc) for E < and Figure (2c) for 
E > 0. From eqns (2.33) and (3.19) to (3.21), the asymptotic form of the density 
profile density is given by 

' 1 + 2E 2 



(3.30) 



_£> 3 (2£) 2 ' 9D 2 

for D » [(1 + 2E)/E 3 ] 1 / 2 and D « [(1 + 2E)/\E\ 3 } 1 / 2 [cf. eqns (3.8) and 
(3.13).] The E > solutions are everywhere underdense relative to the E = 
solutions, going to as (ln|z|) _1 at the origin, whereas the E < solutions are 
everywhere overdense and have a second density singularity. There is a uniform 
core region only in the E = case, although CC show that this also applies for 
E < if there is pressure. 

The form of M(z) = T/S(z) can be deduced immediately from Figure (3a). 
In the E > case it rises or falls monotonically, as in the D = case, so there is 
just one value of z for which M = 1/2. As \z\ — > oo, M tends to a limiting value 

—TO (3 - 31) 

and the apparent horizon will be in z > or z < according to whether this is 
greater or less than 1/2. In the E < case, M will have a minimum where S has 
a maximum. This occurs at 



7rVl + 2£ ' 



(3.32) 



{-2Efl 2 

and since the minimum value is \E\, this is necessarily less than 1/2. One therefore 
has at least two points where M = 1/2, one of which is the apparent horizon for 
the black hole associated with the singularity at zs- Therefore, as in the D = 
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case, there will always be a black hole apparent horizon but not necessarily an 
event horizon. This emphasizes an important difference between the collapse 
singularities at z = 1/D and z = z s : only the latter is associated with an 
apparent horizon, which is why only the former can be naked. 

The D < solutions have the same form as in the D = case, except that 
the E < ones have a second singularity at the value zs given by eqn (3.15) 
with D < 0. In this case, as z goes from —1/D to zs, S first increases to some 
maximum value and then decreases, while V monotonically increases from — oo 
to +oo. As in the E — case, such models are probably physically unrealistic 
since the density diverges due to shell-crossing. They are therefore not shown 
explicitly. The form of the (unphysical) negative-mass solutions, which only exist 
for E > 0, is indicated by the dotted curve in Figures (3). This is similar to the 
D = case shown in Figures (2) except that asymptotically eqn (3.19) applies 
rather than eqn (3.8). 

4 Conclusion 

We may briefly summarize the results of our analysis as follows. (1) There are 
two families of spherically symmetric self-similar dust models: asymptotically 
flat Friedmann solutions and what we have termed asymptotically quasi-static 
solutions. These all represent inhomogeneous cosmological models in which the 
energy function E is constant. They either expand from a Big Bang or collapse to 
a Big Crunch but the singularity is only at t — for the asymptotically Friedmann 
family. (2) Some of the asymptotically Friedmann models represent overdensities 
in a Friedmann background which recollapse to a second singularity and contain 
a black hole which grows as fast as the Universe. The black hole always has 
an apparent horizon but not necessarily an event horizon. Other asymptotically 
Friedmann models represent underdensities in a Friedmann background which 
grow as fast as the Universe. (3) The asymptotically quasi-static models can be 
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interpreted as representing either inhomogenenous cosmological solutions (with 
one or two singularities) if one allows both signs of z or self-similar collapse from 
an initially isothermal distibution if one allows just one sign of z, with a uniform 
density core in one regime and a central black hole or naked singularity in the 
other. (4) We have emphasized the relationship between the z > and z < 
solutions. Any particular asymptotically Friedmann solution is confined to one 
sign of z but any asymptotically quasi-static solution necessarily spans both signs. 

In an accompanying paper (Carr & Coley 1999a), it is shown that the spheri- 
cally symmetric self-similar solutions with pressure share many of the qualitative 
features of the dust ones, especially in the supersonic regime. In particular, all 
of the properties (1) to (4) above still pertain. However, it should be emphasized 
that new types of solution arise when there is pressure. For example, there is an 
exact static solution and an exact Kantowski-Sachs solution, as well as families of 
solutions asymptotic to these. There are also asymptotically Minkowski solutions 
for a > 1/5, some of which asymptote to a finite value of z. The inclusion of 
pressure obviously introduces qualitatively new features in the subsonic regime, 
in particular the possible presence of a sonic point. 

In claiming that our classification is "complete" , it should be emphasized that 
our considerations have been confined to similarity solutions of the simplest kind 
(i.e. homothetic solutions in which the similarity variable is z = r/t). However, 
it should be noted that this is not the only type of similarity. For example, 
Carter & Henriksen (1989) have generalized the concept to include what they 
term "kinematic" self-similarity. In this context the similarity variable is of the 
form z = r/t a for a ^ 1 and the solution may contain some dimensional constant. 
Ponce de Leon (1993) has also introduced the closely related notion of "partial 
homothety". It is not yet clear how easily the analysis of this paper can be 
extended to these cases. Finally it should also be emphasized that we have only 
been studying solutions which are homothetic everywhere. The sort of models 
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considered by Tomita (1997), in which one patches a self-similar transition region 
between other non-self-similar regions, is clearly not covered here. 

Acknowledgments 

The author thanks Alan Coley for useful discussions and is grateful to the 
Yukawa Institute for Theoretical Physics at Kyoto University and the Department 
of Mathematics and Statistics at Dalhousie University for hospitality received 
during this work. 



26 



References 

H. Bondi, 1947, MNRAS 107, 410. 

W. B. Bonnor, 1956, J. Astrophys. 39, 143. 

A. H. Cahill and M. E. Taub, 1971, Comm. Math. Phys. 21, 1. 

B. J. Carr and A. A. Coley, 1999a, Phys. Rev. D., in press. 

B. J. Carr and A. A. Coley, 1999b, Class. Quant. Grav., 16, R31. 
B. J. Carr and S. W. Hawking, 1974, MNRAS 168, 399. 
B. J. Carr and A. Yahil, 1990, Ap. J. 360, 330. 
B. J. Carr and A. Whinnett, 1999, preprint. 

B. Carter and R. N. Henriksen, 1989, Ann. Phys. Supp. 14, 47. 
D. Christodoulou, 1984, Commun. Math. Phys. 93, 171. 

I. H. Dwivedi and P. S. Joshi, 1997, Class. Quant. Grav. 47, 5357. 

C. C. Dyer, 1979, MNRAS 189, 189. 

D. M. Eardley and L. Smarr, 1979, Phys. Rev. D. 19, 2239. 
V. T. Gurovich, 1967, Sov. Phys.-Dokl. 11, 569. 

S. W. Hawking and G. F. R. Ellis, 1973, The Large-Scale Structure of 
Spacetime, Cambridge University Press. 

R. N. Henriksen and K. Patel, 1991, Phys. Red. D. 42, 1068. 

P. S. Joshi and I. H. Dwivedi, 1993, Phys. Rev. D. 47, 5357. 

P. S. Joshi and T. P. Singh, 1995, Phys. Rev. D. 51, 6778. 

27 



R. Kantowski and R. Sachs, 1966, J. Math. Phys. 7, 443. 

A. Krasinzki, 1997, Physics in an Inhomogeneous Universe, Cambridge Uni- 
versity Press. 

K. Lake, 1992, Phys. Rev. Lett. 68, 3129. 

R. Larson, 1969, MNRAS 145, 271. 

S. D. Maharaj, 1988, J. Math. Phys. 29, 1443. 

M. V. Penston, 1969, MNRAS 144, 449. 

A. Ori and T. Piran, 1990, Phys. Rev. D. 42, 1068. 

J. Ponce de Leon, 1993, Gen. Rel. Grav. 25, 865. 

A. M. Sintes, 1996, PhD thesis (University of Balearic Islands). 

R. C. Tolman, 1934, Proc. Nat. Acad. Sci. USA 20, 169. 

K. Tomita, 1995, Ap. J. 451, 1. 

K. Tomita, 1997a, Phys. Rev. D. 56, 3341. 

K. Tomita, 1997b, Gen. Rel. Grav. 13, 625. 

Z. C. Wu, 1981, Gen. Rel. Grav. 13, 625. 

Figures 

FIGURE (1). This shows the form of (a) the scale factor S(z), (b) the ve- 
locity function V(z) and (c) the density function jj,t 2 {z) for the asymptotically 
Friedmann dust models, the arrows indicating the direction of increasing time. 
These are described by a single parameter E where E = in the exact Friedmann 
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case: the z > solutions are overdense and collapse to black holes for E < 
(with an event horizon for E > since V m i n < 1) but they are underdense and 
expand forever for E > 0. The z < solutions are just the time reverse of these. 
The dotted curve corresponds to a solution with negative mass; it is probably 
unphysical but relates to the Kantowski-Sachs solution which arises when there 
is pressure. 

FIGURE (2). This shows the form of (a) the scale factor S(z), (b) the velocity 
function V(z) and (c) the density function [it 2 (z) for dust models with E = 
0. Two different values of D are shown in (a) and (b) but only one in (c). 
These solutions necessarily span both positive and negative values of z. For 
D > they represent monotonically expanding (solid) or collapsing (broken) 
solutions and the latter contain a naked singularity (V m i n < 1) if D exceeds some 
value D + (as assumed here). The D < models (dotted) undergo shell-crossing 
before encountering the singularity and are probably unphysical since V and \x 
go negative. 

FIGURE (3). This shows the form of (a) the scale factor S(z), (b) the velocity 
function V(z) and (c) the density function jit 2 {z) for the asymptotically quasi- 
static dust solutions. These are described by two parameters (D and E) but we 
assume that D is fixed and not all the solutions in (a) and (b) are shown in (c). 
For E > the solutions resemble those in Figure (2), with both monotonically 
expanding (solid) and collapsing (broken) solutions. The collapse singularity is 
naked (V min < 1) if E is less than a value E + (D). For E < there are also 
solutions which recollapse to a black hole (solid) or emerge from a white hole 
(broken), as in the asymptotically Friedmann case. As E decreases, the last 
solution is the symmetrical one for which E = E sym and zg — 1/D, so that the 
solid and broken curves coincide. The curves labelled E = + and E = 0~ show 
the qualitative transition as E passes through 0. The dotted curves correspond 
to an (unphysical) negative mass solution. 
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FIGURE (4). This shows the permitted regime for the parameters E and 
D. The curve labelled E sym indicates the symmetric solution and all physical 
solutions must lie above this. The upper broken line gives the transition between 
different asymptotic forms for and V^. The collapsing solutions have a naked 
singularity in the vertically shaded region below the line labelled E + and the black 
hole solutions have an event horizon and a particle horizon in the horizontally 
shaded region above the line labelled E*. These lines intersect on the lower 
boundary. 
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